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Abstract 



The transport of Brownian particles in a slit geometry in the presence of 
an arbitrary two-dimensional periodic energy landscape and driven by an 
I external force or convected by a flow field is investigated by means of macro- 

^ transport theory. Analytical expressions for the probability distribution and 

^ the average migration angle of the particles are obtained under the Fick- 

C^ Jackobs approximation. The migration angle is shown to differ from the 

Q orientation angle of the driving field and to strongly depend on the physical 

properties of the suspended species, thus providing the basis for vector chor- 
al matography, in which different species move in different directions and can be 
X continuously fractionated. The potential of microfluidic devices as a platform 
'— ' for partition-induced vector chromatography is demonstrated by considering 
T— H the particular case of a piece- wise constant, periodic potential that, in equi- 
r>" librium, induces the spontaneous partition of different species into high and 
^ low concentration stripes, and which can be easily fabricated by patterning 
CN physically or chemically one of the surfaces of a channel. The feasibility to 
; separate different particles of the same and different size is shown for systems 
^ in which partition is induced via Ig-gravity and Van der Waals interactions 
O in physically and chemically patterned channels, respectively. 
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1. Introduction 

Separation of the different constituents of a complex sample has long 
been of paramount relevance in many fields of engineering and science. The 
ubiquity of particle separation, for example, has led to the development of 
a vast number of different techniques [1]. The manipulation of chemical 
and biological species at the micro and nanoscale has also received special 
attention in the ongoing effort toward process miniaturization, with nearly 
every advance in our understanding of the effects that are dominant at these 
scales leading to the innovation of different separation devices as surveyed 
in recent reviews [21 [3l E]. An ideal technique would be able to discriminate 
particles based on small differences on any of a wide range of physicochemical 
properties, thus allowing the fractionation of intricate mixtures. In addition, 
continuous operation is generally preferred over a batch process because of 
the ease of operation and higher yields. To facilitate operation and portability 
it is also desirable that the device be autonomous, thus eliminating the need 
for external components. In contrast, most techniques discriminate particles 
based on a single physicochemical property, are inherently batch-processes, 
or require external fields to drive separative displacement. 

One of the most versatile separation methods is Field-flow fractionation 
— FFF— p], a family of techniques in which a variety of force fields can be 
used to induce separation based on different physicochemical properties, with 
subtechniques including gravitational, sedimentation, flow, thermal, and elec- 
trical FFF [HI [3. In FFF, particles are typically transported by a parabolic 
flow between two parallel plates and, at the same time, externally forced in 
the direction perpendicular to the flow toward one of the walls. The velocity 
of the particles is thus governed by their equilibrium particle-wall separation, 
with particles in closer proximity to the wall moving slower, which results in 
the separation of the sample into bands of particles moving in the direction 
of the flow. Dielectrophorectic methods [HI E] are also versatile since differ- 
ences in a broad spectrum of properties result in different dielectric signatures 
allowing, similarly to FFF, the separation of multifarious mixtures. 

Despite the breadth of FFF techniques, they share the common draw- 
back of being fundamentally a batch process. A technique as wide-ranging 



as FFF that allows for continuous operation is Split-Flow Lateral- Transport 
Thin — SPLITT— fractionation [TO]. In this technique, particles are streamed 
near one of the walls of the SPLITT channel and, at the same time, an ex- 
ternal force transports particles selectively in the direction perpendicular to 
the carrier flow into different streams which can then be continuously col- 
lected. On the other hand, SPLITT only allows bi-modal separation in a 
single stage. Other methods also allowing continuous operation have been 
developed by acoustic and optical means. Acoustic forces have been used 
to separate particles according to their density and compressibility into two 
streams at the pressure node and antinode created by a standing wave be- 
tween two parallel walls [TT]. Optical forces have allowed the sorting of nano 
and microparticles of different materials and sizes, as reviewed by Jonas and 
Zemanek [12]. These continuous methods, however, require the integration 
with external components. 

Continuous separation methods that do not require the presence of ex- 
ternal fields include pinched flow fractionation — PFF— [131 HU [TJ], hydro- 
dynamic filtration — HF— p!6], Hydrophoresis [171 [THl [I9], and inertial fo- 
cusing [201 1211 122] . In PFF particles exit a pinched flow region in different 
stream lines when they are initially aligned to one of the walls of the con- 
striction. Similarly, in HF particles are first aligned along the walls of a main 
channel and then sequentially collected according to their size by controlling 
the flow rate in side channels. Hydrophoresis, on the other hand, does not 
require precise flow focusing of the particles. Hydrophoretic separation is 
achieved using an array of obstacles to induce a pressure field responsible for 
the selective displacement of particles of different size. More recently, iner- 
tial lift forces and Dean flows in curved conduits have been used to stream 
particles in precise locations within a flow channel. However, hydrophoresis 
and inertial focusing methods are difficult to implement for the simultaneous 
collection of particles in polydispersed samples. 

Two dimensional —2D— separation is another approach that allows for 
continuous fractionation. The spatial resolution required for continuous op- 
eration in 2D methods is achieved by a combination of transport in one direc- 
tion with selective displacement in the perpendicular direction and results in 
greater resolving power compared to ID techniques [23] . These methods fall 
in the category vector chromatography [2^ [25]— VC— tecniques, where the 
fractionation relies on differences in the average direction in which species 
being separated move. An extension of FFF, for example, combines radial 
and tangential carrier flows between two parallel disks to allow for the contin- 



uous collection of separated bands in different locations of the fractionation 
cell [26]. Other examples of 2D fractionation use a flow field and an exter- 
nal force applied in the direction perpendicular to the flow as the selective 
displacement. In magnetophoresis [27j, for instance, particles are sorted ac- 
cording their magnetic susceptibility and size. 

A promising trend in VC exploits the interactions between the species 
being separated and features intentionally designed in the separation media 
to drive separative displacement. For example, deterministic lateral displace- 
ment [281 ISHl EQI |3T1 [32], [33] — DLD— is carried out in devices patterned with 
a 2D sieving matrix in which particles of different sizes move at different an- 
gles, thus enabling the simultaneous fractionation and collection of multiple 
components in a polydispersed sample with very high resolution. However, 
the underlying mechanism leading to separation in these sieving devices is 
not completely clear, which has led so far to ad hoc designs that are difficult 
to optimize [311 [351 [36] • Iii this context, Dorfman and Brenner [2^ considered 
the illustrative case of a periodic system consisting of repeating layers of two 
inmiscible fluids. They showed that differences in the partition ration of the 
particles —between the two phases—, due to differences in affinity between 
the particles and each of the layers of fluid, results in vector chromatography 
when the particles are animated by a constant external force. Unfortunately, 
as noted by the authors [21], such a two- fluid layer system cannot be easily 
implemented in practice. 

In this study we demonstrate the potential of planar microfluidic devices 
as a platform to achieve partition-induced vector chromatography — PIVC— , 
a technique that combines the versatility of FFF with the multiple advan- 
tages of continuous 2D vector separation. The high surface-to-volume ratio 
characteristic of microfluidic devices makes it possible to induce partition by 
means of surface interactions. In fact, it has been shown that the sponta- 
neous partition of different species can be controlled by means of the energy 
landscape created by a chemical or physical pattern on one of the surfaces 
of a channel [371 EH]- In particular, we consider cases in which partition is 
induced via Van der Waals forces and 1-g gravity, both approaches resulting 
in autonomous devices. In addition, external flelds could also be used to 
cause or enhance partition. We show that partition results in diffusive fluxes 
that are responsible for the migration of particles at angles different from the 
orientation angle of the driving fleld. 

The content of this work is structured as follows. First, the equation 
governing the probability density is derived for a spatially periodic potential 



following the macrotransport paradigm [3S] • In the subsequent sections, this 
governing equation is solved under the Fick- Jacobs — FJ— approximation 
for the cases in which the particle is animated either by a constant external 
force or by a fluid flow. Closed-form solutions for the particle trajectory angle 
in a slit geometry are obtained for an arbitrary two-dimensional potential. 
These general solutions are evaluated to obtain analytical expressions for 
a piece-wise constant potential. Then, it is shown how such a potential 
can be achieved by patterning chemically or physically one of the surfaces 
of a microfluidic device with an array of rectangular stripes. Specifically, 
we consider the case in which partition between the stripes is induced by 
differences in the potential energy resulting from 1-g gravity. Van der Waals, 
and electrostatic forces. Finally, we look at experimentally accessible systems 
and validate the results obtained with the FJ approximation with Brownian 
dynamics simulations to show that particles exhibiting different partition 
ratios can be effectively fractionated via vector chromatography. 



2. Particle transport in patterned niicrofiuidic devices: 
transport theory 



Macro- 



In this section we derive, by means of macrotransport theory [39], the 
macroscopic equations governing the transport of suspended particles in the 
periodic system shown in Fig. [l] The particles are confined between two 
infinite parallel walls separated by a distance d. We shall assume that the 
physicochemical properties of the walls of the channel render the system Ix- 
and /y-periodic in the x and y directions, respectively. Let P(i2, t|i2o) be 




Figure 1: The particle is suspended between two parallel walls separated a distance d —slit 
geometry— and the physicochemical properties of the impermeable walls are such that the 
system is periodic in the x and y directions. 



the conditional probability density of finding the particle center, modeled as 



a stochastic variable, at location R = (x, y, z) at time t, given the initial 
position Rq = (xq, 70,2^0) at time t = 0. This conditional probability density 
is governed by the following convection diffusion equation, 

dP 

— + V-J=6{R-Ro)6{t), (1) 

where the right-hand side represents an instantaneous unit impulse at Rq 
at time t = 0. The flux J is composed of a convective and a diffusive 
contribution, the latter assumed to have Fickian form, 

J = UP~D-VP, (2) 

where U is the net convective velocity of the particle due to both flow and 
applied forces, and D is the diffusion tensor. Due to the periodicity of 
the system, P(i2, t|i2o) can be written in the equivalent functional form 
P{Rn^m, ■*", ^|-Rno,mo) ^o) wherc Rj^k = j Ix ^x + k ly iy represents the location 
of the {j, k)th unit cell with respect to some arbitrary origin, and r = [x, y, z) 
is the intracell position vector. In addition, both U{r) and D{r) are periodic 
and dependent on the intracell position r only. Therefore, the probability 
distribution function P does not depend on Rn,m and Rno,mo independently 
but on their difference, P{Rn,m — Rno,mo,'f,'t\^o)- Then, the zeroth order 
moment of the distribution, 

Po{r,t\ro) = ^P{Rn^rn - Rno,mo, r,t\ro), (3) 



represents the conditional probabihty density that the particle will be found 
in the local, intracellular position r at time t, given that it was introduced 
at the location Tq at time t = 0. PQ{r,t\rQ) satisfies the equation 

BP 

-^ + V-Jo = 5{r-ro)6{t), (4) 

where the flux Jo follows Eq. (|2]). This governing equation for the proba- 
bility distribution corresponds to the Smoluchosky equation describing the 
stochastic motion of a particle in the limit of negligible inertia |10]. The 

steady state solution, 

Po°^(r)=limP(r,t|ro), (5) 

I— ^-oo 



is then governed by the equation 

V ■ J^ = 0, (6) 

with periodic boundary conditions, 

P^{x,y,z)=P^{x + h,y,z), (7a) 

P^{x,y,z) = P^{x,y + ly,z), (7b) 

and the no-flux condition enforced on the impermeable walls confining the 
particle in the z-direction. Finally, the particle is bound to be somewhere 
within the volume r of the unit cell, which implies the standard normalization 
condition 

'' P^{r)d^r = l. (8) 



The asymptotic flux of probability density can be written in our case as 

J^ = uP^ + M- FP^ - M- VVP^ - D ■ VPo°°, (9) 

where u{r) is the velocity of the particle when convected by a flow field, F is 
a constant external force acting on the particle, V{r) is the potential energy 
landscape in which the particle is moving, and ]V[{r) is the mobility tensor, 
and all quantities are assumed to be periodic. For the geometry considered 
here, it is convenient to write the mobility and diffusion tensors as 

M=44M^ + (/-44)M||, (10a) 

D=iAD^ + {I-iX)D\\. (10b) 

The mobilities and diffusivities follow the Stokes-Einstein relation, D^ = 
ksTMx and D\\ = kBTM\\, where ks is the Boltzmann constant and T the 
absolute temperature. 

Knowledge of P^^ yields the mean particle velocity from the quadrature 

If = [ J^dh. (11) 



The orientation angle of if is the relevant parameter in vector chromatog- 
raphy, termed hereafter the chromatrographic trajectory angle, more simply 
the trajectory angle |23], or the migration angle. In the following sections 
we derive the trajectory angle under the Fick- Jacobs approximation for the 
cases in which the particle is driven either by an externally applied force or 
by a flow field. 



2.1. Non-equilibrium transport: The Fick- Jacobs approximation 

The three-dimensional problem posed in the preceding section can be 
solved numerically for a particular geometry and driving fields. In this sec- 
tion, we focus on essentially 2D problems that capture the fundamental mech- 
anisms leading to PIVC, and which are amenable analytically. Specifically, 
we consider the x-invariant case whereby u{r), M{r), V{r) and Pl^{r) de- 
pend on y and z only, and the periodicity in the x-direction is automatically 
satisfied. In addition, we study 2D driving fields in which the external force 
and the fiuid-fiow velocity far from the particle are in the x-y plane. Under 
these conditions the asymptotic fiux, Eq. ^, takes the form 

[ dV dP°^ \ 

+ i^UyP,- + M\\FyP^ - M\\—P^ - ksTM\\ ^ J iy (12) 



+ ( -mJ^P^ - knTM ^^' 



dz '' "^ ^ dz 

and Eq. ([6]), governing the asymptotic probability density, becomes 



d / dV dP°^ 

A heuristic approximation first introduced by Jacobs [HI |12] allows solv- 



ing Eq. (13) analytically assuming local equilibrium in the transverse di- 
rection. This approximation, termed by Zwanzig the Fick- Jacobs (FJ) ap- 
proximation, has been shown to provide accurate results in narrow chan- 
nels [l3l HU HH]. Under this local equilibrium assumption, and given that 
the system is a;-invariant, the asymptotic probability distribution conditioned 
to a given intracelular position y can be written as 

pix, z\y) ^ p{x, z\y%^ = e-^MI^^T j j^y^^ (^4) 

where 

I{y)= [ [ e-^^y''^/''^^dzdx. (15) 



Therefore, the position of the particle in y suffices to determine the distribu- 
tion, which takes the form 

Po^ir) = V{y)p{x,z\y) ^ V{y)p{x, z\y),, = V{y)e'''^''^^/'-^ / I{y), (16) 



where 



p(.)^// 



P^{r)dzdx 



(17) 



is the marginal probabihty density. The cases in which the particle motion 
is driven either by a constant external force or by a fluid flow are considered 
under the FJ approximation in what follows. 

2.1.1. External force 



When the particle is animated solely by a constant external force, Eq. ( 13 ) 
reads 

dy 

+ ^ ( -M^^P^ - ksTM^^^^ ) . (18) 







M\\FyP^ 



Mil 



-P^ - ksTMi 



dz 



dV 



Qpoo 

dP^' 



dz 



Nondimensionalizing with the variable changes x — >■ x/lx, y — )• y/ly, z — )■ z/d, 
M|| -^ M||/Moo and Mj_ -^ M^/M^o -where M^o = 67r/ia being the mobihty 
of an unbounded particle of radius a in a fluid of viscosity /i— , V — )■ V/kBT, 
and P^ — !■ tP^, substituting the FJ approximation for P^, integrating in 
X and z, and rearranging the terms yields 

d 







dyj 


J J^^{r)dzdx 








df 
dy\ 


1 1 M\\{y,z)e-^^y^'^dzdx 


(Pe^(^)- 

l m 


d 

dy 


\V{y)] 
[i{y)\ 



(19) 



where the only dimensionless parameter is the Peclet number, Pe = Fyly/ksT. 
The total dimensionless flux in the y-direction, J, is clearly constant and can 
be written as 



J 




J^^{r)dzdx = I{y)M\\{y) 



VeV{y) 



dy 



where 



my) 




M\\{y,z)e~^^y''^dzdx / liy), 



(20) 



(21) 



is the local average of the mobility, and Vi^y) = V{y)/ I{y). The general 
solution for the marginal density is given by 

dy' 



V{y) 



,Pej/ 



-J 



e^-y'I{y')M\\{y'] 
9 



+ N 



m- 



(22) 



The flux J and the constant A^ are obtained by solving simultaneously the 
periodicity and normalization conditions: 



A^ 



dy' 



1 



1 n+i 

dye'^'yiiv) / 

Jy 



-Po 



) 



dy' 



e^-y'I{y')M\\{y'y 
The total flux in the z direction is identically equal to zero. Thus, 



If = u: 4 + u; iy 



Ox '^ ^ '^x 



'^Qy d r ly 



(23a) 

(23b) 
(23c) 

(24) 



Employing Eqs. (23) and (24), we obtain one of our main results, a chro- 



matographic trajectory angle that depends in general on the properties of 
the particle through the energy landscape V and the external force F: 



U* Pe 

tan e* = ^ = tan Op 

U* 
y 



1 — p-Pe , 



dye'''yi{y)M\\{y) 



y+l 



dy' 



,y I{y')M\\{y')e^-y" 
(25) 
where 9p = aictajLi{Fx/Fy) is the orientation angle of the external force. 



2.1.2. Entrainment in a fluid flow 



When the particle is convected by a fluid flow, Eq. (13) reduces to 
d 







dV dP°° 



dz 



oz oz 



(26) 



Nondimensionalizing with the same change of variables as those used for 
Eq. (18) in addition to Uy — )■ Uy/U, where U is the average velocity of the 



particle- free flow, substituting the FJ approximation for P^, integrating in 
X and z, and rearranging the resulting terms we obtain the —constant— total 
flux in the y direction 



J 




J^{r)dzdx = I{y)M\\{y) 



Pg ^y^y) ^ 



M\\{y) 



Hy) 



dvjy) 
dy 



(27) 
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In this case, the Peclet number is Pe = Uly/kBTMoo and Uy{y), analogously 
to M\\{y), is the local average of the velocity Uy. We can then define a local 
effective force as the ratio of the local average velocity to the local average 
mobility, 



<P{y) 



M\\{y)' 



(28) 



Then, the general solution for the marginal probability distribution under 
the FJ approximation can be expressed as 



nv) 



,Pe $(j/) 



-J 



dy' 



eP^'^(y')Iiy')M\\{y': 



+ N 



-V{y,z) 



where 



^{y) 



(j){y')dy'. 



(29) 



(30) 



The flux J and the constant A^ are obtained, as in the external force case, by 
solving simultaneously the periodicity and normalization conditions: 



N 



dy' 



1 



J 



[1 



-Pe$(l) 



dye^'^^'^y^Iiy) 



J/+1 



dy' 



3Pe<&(y')/(^')M||(?/0 



Lastly, the chromatographic trajectory angle is given by 



tanr 



U* 

Tl* 



tan 6* 



Pe 



l_e-Pc*(i) 



rf2/e^^*^^^J(2/)M||(y) 



y+l 



(31a) 

(31b) 

(31c) 

dy' 



y eP-'^(y')I{y')M\\{y') 
(32) 



where 9f is the orientation angle of the particle- free flow. Eqs. (25) and (32) 
are valid for an arbitrary two dimensional potential, under the FJ approxi- 
mation. These expressions are used in the next section to obtain analytical 
results for a potential induced by surface modiflcations in a microfluidic de- 
vice. 
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3. Vector chromatography in microfluidic devices: Stripe patterns 

Using standard microfabrication techniques, a straight channel in a mi- 
crofluidic device can be modified to yield a periodic system by patterning one 
of the surfaces of the channel with an array of rectangular stripes as shown 
in Fig. [2J We consider the case where the interfaces between the stripes are 
sharp, resulting for particles at close proximity to the wall, in approximately 
a piece-wise potential of the form 

Vi{z) 0<y<e 
V2{z) e<y <l. 



V{r) 



(33) 
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Figure 2: Schematic view from the top of a particle driven by an external force. The 
particle is suspended over a pattern of rectangular stripes of widths e and 1 — e relative to 
the period length ly. 

The migration angles can then be calculated by substituting this potential 



into the general expressions given by Eq. (25) and (32). It is worth mention- 



ing that the transition potentials do not contribute to the trajectory angle 
whenever the transition regions are much narrower than the stripes, and 



therefore their explicit inclusion in Eq. (33) is not necessary. This simple ap 



proximation provides analytical expressions that facilitate the understanding 
of the physics underlying the separation process, and shows that separation 
is feasible. 
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3.1. Interaction potential 

In general, the potential on each stripe includes the Van der Waals and 
electrostatic forces between the particle and the surface and, for non-neutrally 
buoyant particles, it has a gravitational contribution as well, 

V{h) = Vedi{h) + V^Uh) + Vgrav{h), (34) 

where h is the particle wall separation —see Fig.[T}-. The electrostatic double- 
layer potential for the case of a z-z electrolyte has the form 



with 



Vedi{h) = aBexp{-Kh) (35) 



where e^ is the dielectric permitivity of the medium, e is the charge of an 
electron, ipp and ips are the Stern potentials of the particle and the wall, 
respectively, and 

K~' = f ^eVaiVA/eA;Br J (37) 

is the Debye length, where Cj is the bulk electrolyte concentration of species 
i, z is the valence of the ions, and Na is Avogadro's number. 

The Van der Waals potential can be rigorously computed via Lifshitz 
theory and fit for convenience to a power law expression of the form [37] 

V^dUh) = -aAh-P. (38) 

These surface potentials have been referenced to zero at an infinite separation. 
The gravitational potential is given in terms of the buoyant weight of the 
particle 

Vgravih) = A/STia^Apg {h + Kef)/ kBT, (39) 

where Ap = Pp — Pm is the density difference between the particle and the 
medium —or buoyant density—, g is the acceleration due to gravity, and href 
is the height of channel wall with respect to an arbitrary reference. Fig. [3] 
shows the total potential for particles with a radius of 2 pm as a function of 
the particle-wall separation, calculated using parameters representative of the 

13 



experimental values reported by Wu et al [57] for two values of the buoyant 
density. Specifically, Stern potentials of -60 mV, and Van der Waals afinity 
prefactor and exponent of of 2 nm and 2, respectively, are representative of 
the experimental values corresponding to a silica particle suspended above a 
bare glass surface in a 1 mM aqueous solution of a 1:1 electrolyte. In this 
case the surface interactions are important within a few hundred nanometers 
from the wall with gravity dominating for larger separations. The less dense 
particle is confined to a broader range of separations and is, on the average, 
farther away from the wall than the denser particle —see panels (a) and (b) 
in Fig. |3} respectively. Fig. |4] shows the average particle-wall separation and 
the standard deviation as a function of particle radius for different buoyant 
densities. The particles are on average closer to the wall and more narrowly 
confined as the particle buoyant density and size increase, with the excur- 
sions of the heavier particles restricted within tens of nanometers from their 
equilibrium position. This degree of confinement makes the FJ approxima- 
tion valid and in fact accurate even at relatively large Peclet numbers, as 
shall be shown later by comparison with Brownian dynamic simulations. 
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Figure 3: Total potential in fc^T units —left axis, light shaded region— and corresponding 
Boltzmann distribution —right axis, dark shaded region— as a function of the particle- 
wall separation for (a) a particle twice as dense and (6) slightly denser than the aqueous 
medium. The inset in (a) shows the repulsive energy barrier that has to be overcome to 
attain the primary minimum. The particle has a radius of a = 2 fj,m, the parameters of 
the electrostatic potential are "0^ = 'i/'s = — 60 mV with an aqueous molarity of 1 mM —1:1 
electrolyte—, and the parameters of the Van der Waals potential are A = 2 nm, p — 2. 
The surface and gravitational potentials vanish at an infinite separation and at the wall, 
respectively. 
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Figure 4: Average particle separation —solid curve— and corresponding standard deviation 
—dashed curve—. The standard deviation is comparable to the average separation for the 
particles in the low range of buoyant weights considered. As the buoyant density and 
the size of the particle increase, the particles are closer to the wall and more narrowly 
confined. The parameters of the electrostatic and Van der Waals potentials are as noted 
in the caption to Fig. [3] The arrow points in the direction of increasing buoyant density. 
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3.2. Inducing partition via physical and chemical patterns 

In the absence of a driving field, the concentration of particles undergo- 
ing Brownian motion over a stripe pattern as shown in Fig. [2] follows the 
Boltzmann distribution. In this case, the marginal distribution is given by 

<y < e 



Viy) = < '^^'r^^ _,_, , (40) 



e+(l-e)K 

with the partition ratio, K, defined as 



e<y<l. 



K= e-^'^'^dz e~^'^'^dz. (41) 



It is clear that differences in the potential energy of the particle in the two 
regions would result in partition, K being a measure of the relative affinity 
between the particle and each of the stripes. 

Patterning the stripes chemically would result in different surface po- 
tentials while patterning the stripes physically modifying the topography of 
the surface would result in different sampling of the gravitational potential. 



Fig. |5j shows partition ratios predicted from Eq. (41) with the electrostatic 
and Van der Waals potentials based on the experimental values reported by 
Wu et al [33 • Specifically, we first show the partition ratio induced by grav- 
ity in a pattern of chemically identical stripes as a function of the height 
difference between them for —Fig. |5](a)— different values of the buoyant den- 
sity, ranging approximately from that of polystyrene particles to that of silica 
particles, and —Fig. Islfe)- for different values of the particle radius and for 



a buoyant density representative of a silica particle. In this case, in which 



the surface properties of both stripes are the same, Eq. (41 ) for the partition 
ratio simplifies to 

K = e^/STraM-Hzlpg/fcsr ^2) 

where ATi is the height difference between the stripes. The partition ratio 
increases with the height difference between the stripes, the buoyant density 
and the size of the particle. It is worth mentioning that the sampling of the 
transverse direction is such that the average particle-wall separation is larger 
than the largest AH considered —see Fig. |4]— , and therefore the particles 
do not experience physical hindrance at the transition regions. In panels (c) 
and (d) of Fig. [slwe present partition induced via Van der Waals interactions 



by increasing the particle affinity, parameter A in Eq. (38), with one of two 

16 



—otherwise identical— stripes. Specifically, a value of p of 2 and values of A 
between 2 and 10 nm are representative of the range spanned in the case of 
a silica particle in a 1 mM aqueous solution of a 1:1 electrolyte interacting 
with a glass substrate coated with a gold film of thickness between and 30 
nm. In this case the Van der Waals-induced partition ratio increases as the 
affinity ratio, the buoyant density, or size of the particle increase. 

3.3. Non-equilibrium transport 

We now consider the cases in which the particle is either driven by a 
spatially uniform external force F oriented at an angle 6p with respect to the 
y axis —see Fig. [2]—, or convected by a flow field that far from the particle is 
oriented at an angle 6f with respect to the y axis. The migration of particles 
at angles that differ from the forcing direction can be explained qualitatively 
by considering the normal flux at the transition between the stripes. Consider 
the equilibrium marginal probability density shown in Fig. |6](a) for a partition 
ratio K = 2. In the presence of a driving field, the convective fiux on each 
side of the interfaces between the stripes would differ due to the concentration 
differences. Thus, for the system to reach steady state and the corresponding 
continuity of the total flux, concentration gradients and the ensuing diffusive 
fluxes must compensate for the differences in the convective fluxes. As can 
be seen in Fig. |6] {b)-{d), the diffusive fluxes augment —reduce— the flux in 
the region of lower —higher— concentration and convective transport. The 
average diffusive flux on each stripe is proportional to the concentration 
difference across them —height of the shaded regions shown in Fig. [6]—. We 
also show in Fig. [6] that the average diffusive fluxes differ by a factor K. 
The net result is a diffusive component of the flux normal to the stripes that 
reduces the total flux in this direction, the effect being stronger for increasing 
partition ratios. As a consequence, the trajectory angle of the particle is 
different from the orientation angle of the fleld and strongly depends on the 
partition ratio. This is the basis for the separation of particles exhibiting 
different partition ratios. 

Evaluating the general expression for the chromatographic trajectory an- 



gle, Eq. (25), for the piece- wise potential given by Eq. (33), we obtain the 
trajectory angle in terms of the partition ratio for the case in which the 
particles are driven by an external force 

(l-Air)Ml-e-P-) n -e-P-^^i-^)] 1 , 
tanr = tan^^<il + ^^^^ ^^-L^ i_ ^ , (43) 
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Figure 5: Partition ratio induced —(a) and {b)~ by gravity in a pattern of chemically 
identical stripes as a function of the height difference between the two stripes, and —(c) 
and (d)— by Van der Waals interactions as a function of the particle-surface affinity with 
one of two otherwise identical stripes. In (a) and (c) the partition ratio is plotted for 
different values of the buoyant density, whereas in (6) and (d) the partition ration is 
plotted for different particle sizes. The arrow traverses curves of increasing value of the 
respective parameter. For all cases, the Stern potential of the particle and stripes is -60 
mV, the molarity of the aqueous solution —1:1 electrolyte- is 1 mM , and the Van der 
Walls exponent is p = 2. 
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Figure 6: Marginal probability density for different Peclet numbers for a partition ratio 
K = 2. The differences in the convective fluxes at the interface between the stripes due 
to partition result in compensating diffusive fluxes to satisfy continuity of the flux. The 
shaded areas represent the concentration difference across the stripes, which differ by a 
factor of K. The net diffusive flux opposes the total flux in the direction normal to the 
stripes. 
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where A = M\\^/ M\\^ with M\\. the local average of the mobility over stripe 
i defined by Eq. (21). In the limit of vanishing Pe the above expression 
reduces to 



tan^* = tan^j 



1 + 



\Ky 



\K 



efl 



(44) 



while in the convection dominated case the trajectory angle tends to the 
forcing angle. 



Analogously, replacing the piece- wise potential into Eq. (32) yields the 



migration angle for the case in which the particle is convected by a fiuid fiow 

:i - XK) (1 - nXK) (1 - e-P^-^i^) [1 - e-P^<^in(i-^)] 1 



tan^* = tajidfi 1 



Ii\K 



-Pc<^l[e+n(l-e)] 



Pe0i 
(45) 
where 11 = 02/</'i; with (/)i = Ui/ M\\. being the effective hydrodynamic force 
acting on the particles moving on stripe i. 

In the limit as the Pe tends to zero the above expression simplifies to 



tan 6* = tan 9f 



'l-\K){l-U\K) e(l-e)n 



XK 



nil 



(46) 



and as in the external force case, the migration angle converges to the forcing 
angle in the high Pe number limit. 

3.4- Brownian dynamic simulations 

The expressions for the migration angle in the case of stripe patterns, 
Eqs. (43) and (45), are valid under the FJ approximation, i.e., under the 
assumption of instantaneous equilibrium in the transverse direction. This 
assumption is clearly valid when the particle has ample time to equilibrate 
in the direction perpendicular to the patterned surface via diffusion as it is 
convected across the stripes by either an external force or by a fluid flow. 
Therefore, the relevant parameter is the ratio between the diffusive time 
in the transverse direction and the convective time along a single cell, i.e., 
Pe = Pe {(j/ly) where a is the width of the conflnement region. Then, 
if the particles are narrowly conflned to the vicinity of the patterned sur- 
face the FJ approximation could be valid even for large values of the Peclet 
number P3| H5]. The particles could be physically constrained to move in a 
narrow channel or, alternatively, a highly conflning potential could restrict 
the vertical motion of the particles to a narrow region independently of the 
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height of the channel. Therefore, in order to determine the full range of va- 
lidity of the results presented in the preceding section and, more important, 
the range of Pe that result in separation, Brownian dynamic — BD— simula- 
tions were carried out for the case when the particle is convected by a fluid 
flow. Using the Ermak and McCammon algorithm |47], the dimensionless 
equations of motion for a suspended particle are 



Ax = CPe v° sin Of At + ^J2Ajj Z\t x (47a) 

fdvV 



Ay = Pe v^ cos Of At -^11(7^) ^t + J2AjjZ\t ip (47b) 

Az = 7^ (^] ° At - 7'A^ (^\ ° At + 7^2A^Z\t co (47c) 

The characteristic scales used to nondimensionalize the equations are l^, ly, 
and d in the corresponding directions, the unbounded mobility —Moo— and 
diffusivity —Doc— for the corresponding quantities, the average velocity of 
the particle- free flow —U— for the particle velocity, fc^T for the interaction 
potential, and ly/Doo for the time step. In the above equations ^ = Ix/ly, 7 = 
Iz/ly, V is the dimensionless particle velocity —in the forcing direction— due 
to the flow, A is the dimensionless mobility —or diffusity— , Pe = Uly/Dao, 
Of is the orientation angle of the particle free flow, and Xi '4'i and u are 
independent random numbers chosen from a symmetric distribution of zero 
mean and unit width. Lastly, the superscript indicates that the quantities 
are evaluated at the beginning of the time step. 

In order to carry out the simulations it is necessary to specify transition 
potentials. As the particle traverses the transition region from one stripe to 
the other, it experiences a force toward the stripe of higher affinity. The mag- 
nitude of this force along the transition region is dictated by the imbalance 
in the particle-substrate affinity between the anterior and posterior regions 
of the particle. The exact form of the transition potential would depend on 
physico-chemical properties of the transition and these in turn would depend 
on the specific way in which the pattern was created. Exact expressions for 
electrostatic and Van der Waals interaction potentials are not available even 
for the case in which the properties of a half-space change discontinously at 
a point resulting in two regions of homogeneous but different properties. We 
model the transition potential as an s-shape curve given by a third degree 
polynomial in y constructed to satisfy continuity of the total potential and 
yield a vanishing force at the edges of the transition regions at every height. 
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3. 5. PIVC in microfluidic devices 

To show the potential of PIVC we consider different separation scenarios 
of a mixture of particles carried by a pressure-driven flow through the stripe 
pattern shown in Fig. |2] and compare the predictions based on the FJ ap- 
proximation with the results of BD simulations. The geometry considered is 
a slit channel with a height of 10 /im. The bottom wall is patterned with 
10 /im stripes, oriented at a forcing angle of 45 degrees with respect to the 
average flow —6f = 45°—. To predict the migration angles from Eq. (45) we 



use the analytical approximation for the particle mobility given by Pawar 
and Anderson |48j and interpolate the numerical results tabulated by Staben 
et al. [l9] for the velocity of a particle confined between two parallel walls. 
Let us note that in the case of physical patterns we neglect the effect of the 
steps in the hydrodynamic force and resulting velocity of the particle. Fig. [7] 
presents the discrimination angle, i.e., the difference between the trajectory 
angle and the forcing angle, for different pairs of particles as a function of 
the average velocity of the carrying flow. Specifically, each of the systems 
and pair of particles considered in panels {a)-{d) in Fig. [7] correspond to a 
system and particles for which the partition ratio is shown in panels {a)-{d) 
in Fig. [5} respectively. The correspondence is indicated by solid squares and 
stars. In Fig.JTJ^a) —{b)— particles of the same —different— size are separated 
according to their mass with partition induced by gravity. Panels (c) and 
(d) in Fig. [T] show the separation of particles based on their Van der Waals 
interaction with chemically different stripes. Fig. W[c) shows the separation 
of particles of the same size but with different affinity with one of the stripes. 
The separation of particles of the same material but of different sizes is shown 
in Fig. ^d). Separation is readily accomplished for all cases considered; one 
of the particles only marginally departs from the forcing angle owing to weak 
partition, while the other particle is restrained to one of the stripes resulting 
in strong deviations from the forcing angle. The data points and respective 
error bars correspond to the average and standard deviation, respectively, of 
128 realizations of BD simulations tracking the particle for 50 units of di- 
mensionless time and with a transition region with a width of two-hundredth 
the length of the unit cell in the ^/-direction. The FJ approximation and the 
simulation results are in excellent agreement over the whole range of discrim- 
ination angles due to the narrow confinement of the particles by the potential 
—see Fig. |4]-. 
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Figure 7: Discrimination angles as a function of the average velocity of the underlying 
flow with partition induced by gravity, (a) and (6), and Van der Waals interactions, (c) 
and (d), for particles of the same and different sizes. The corresponding partition ratios 
are indicated by the solid square and star in the plots in Fig. [5] with the same labels. 
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4. Conclusions 

We showed that vector chromatography in planar microfiuidic devices 
is feasible by harnessing surface interactions that lead to the spontaneous 
partition of different species. First, we derived analytical solutions for the 
trajectory angle of a particle driven by an external force or convected by a 
flow field in a slit geometry for an arbitrary two-dimensional periodic po- 
tential under the Fick- Jacobs approximation. We showed that the migra- 
tion angle depends on particle properties, thus providing the basis for vector 
chromatography of different species. We considered the case of a piece-wise 
constant periodic potential that could be created by chemically or physi- 
cally patterning one of the surfaces of a microfiuidic device with an array 
of rectangular stripes, thus causing the spontaneous partition of suspended 
particles. We showed that partition results in diffusive fluxes that oppose the 
total flux and reduce the velocity component in the direction normal to the 
stripes making the particles migrate at angles that differ from the orientation 
angle of the driving field, the effect being stronger with increasing partition 
ratios. Specifically, we considered the fractionation of particles of the same 
and different size using both physically patterned channels where partition 
is induced by 1-g gravity as well as chemically patterned channels in which 
partition is induced by Van der Waals interactions. We validated our results 
by means of BD simulations that agree well with the results obtained by 
means of the FJ approximations due to the highly confining nature of the 
particle-wall interaction. Thus, we demonstrated that partition-induced vec- 
tor chromatography — PIVC— of particles exhibiting different partition ratios 
is feasible in microfiuidic devices. In addition to continuous operation, an- 
other promising feature of PIVC is its potential versatility that stems from 
the vast range of physicochemical properties that can be employed to in- 
duce partition. Furthermore, separation can be achieved without the need 
for external components making the device autonomous, thus facilitating 
its portability and integration with other lab-on-a-chip components. In ad- 
dition, partition can also be induced or enhanced by means of externally 
applied fields. As in FFF, different subtechniques are envisioned depending 
on the nature of the forces use to cause partition. Further work is needed to 
explore different pattern geometries besides the stripe case considered here in 
order to optimize the sensitivity of the discrimination angle to the partition 
ratio and extend the range of operation to higher Peclet numbers. Lastly, it 
is interesting to note that the results presented here suggest the possibility 
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of inducing vector chromatography by exploiting the dependence of the tra- 
jectory angle on the velocity and mobility of the particles over each stripe 
through the parameters A and 11. 

[1] J. C. Giddings, Unified Separation Science, A Wiley-Interscience publi- 
cation, USA, NJ, 1991. 

[2] N. Pamme, Lab Chip 7 (2007) 1644-1659. 

[3] T. Kulrattanarak, R. Van der Sman, C. Schroen, R. Boom, Adv. Colloid 
Interface Sci. 142 (2008) 53-66. 

[4] M. L. Kovarik, S. C. Jacobson, Lab Chip 81 (2009) 7133-7140. 

[5] J. C. Giddings, Science 260 (1993) 1456-1465. 

[6] J. C. T. Eijkel, A. van der Berg, Lab Chip 6 (2006) 19-23. 

[7] T. Kowalkowski, B. Buszewski, C. Cantado, F. Dondi, Grit. Rev. Anal. 
Chem. 36 (2006) 129-135. 

[8] R Gascoyne, J. Vykoukal, Electrophoresis 23 (2002) 1973-1983. 

[9] M. R Hughes, Electrophoresis 23 (2002) 2569-2532. 

[10] J. C. Giddings, Sep. Sci. Technol. 20 (1985) 749-768. 

[11] M. Evander, A. Lenshof, T. Laurell, J. Nilsson, Anal. Chem. 80 (2008) 
5178-5185. 

[12] A. Jonas, R Zemanek, Electrophoresis 29 (2008) 4813-4851. 

[13] M. Yamada, M. Nakashima, M. Seki, Anal. Chem. 76 (2004) 5465-5471. 

[14] J. Takagi, M. Yamada, M. Yasuda, M. Seki, Lab Chip 5 (2005) 778-784. 

[15] Y. Sai, M. Yamada, M. Yasuda, M. Seki, J. Chromatogr. A 1127 (2006) 
214-220. 

[16] M. Yamada, M. Seki, Lab Chip 5 (2006) 1233-1239. 

[17] S. Choi, J. Park, Lab Chip 7 (2007) 890-897. 

[18] S. Choi, J. Park, Lab Chip 9 (2009) 1962-1965. 

25 



[19] S. Choi, S. Song, C. Choi, J. Park, Anal. Chem. 81 (2009) 50-55. 

[20] D. D. Carlo, D. Irimia, R. G. Tompkins, M. Toner, Proc Natl Acad Sci 
USA 104 (2007) 18892-18897. 

[21] A. A. S. Bhagat, S. S. Kuntaegowdanahalli, I. Papautsky, Lab Chip 8 
(2008) 1906-1914. 

[22] J. Park, S. Song, H. Jung, Lab Chip 9 (2009) 939-948. 

[23] J. C. Giddings, Anal. Chem. 56 (1984) 1258A-1270A. 

[24] K. D. Dorfman, H. Brenner, J. Colloid Interface Sci. 238 (2001) 390-413. 

[25] K. D. Dorfman, H. Brenner, Phys. Rev. E 65 (2002) 052103. 

[26] P. Vastamaki, M. Jussila, M. Riekkola, Analyst 130 (2005) 427-432. 

[27] N. Pamme, A. Manz, Anal. Chem. 76 (2004) 7250-7256. 

[28] L. R. Huang, E. C. Cox, R. H. Austin, J. C. Sturm, Science 304 (2004) 
987-990. 

[29] J. A. Davis, D. W. Inghs, K. J. Morton, D. A. Lawrence, L. R. Huang, 
S. Y. Chou, J. C. Sturm, R. H. Austin, Proc Natl Acad Sci USA 103 
(2006) 14779-14784. 

[30] P. Mao, J. Han, Lab Chip 9 (2009) 586-591. 

[31] Z. Li, G. Drazer, Phys. Rev. E 65 (2007) 050602. 

[32] K. J. Morton, K. Loutherback, D. W. Inghs, O. K. Tsui, J. C. Sturm, 
S. Y. Chou, R. H. Austin, Lab Chip 8 (2008) 1448-1453. 

[33] K. J. Morton, K. Loutherback, D. W. Inghs, O. K. Tsui, J. C. Sturm, 
S. Y. Chou, R. H. Austin, Proc Natl Acad Sci USA 105 (2008) 7434- 
7438. 

[34] J. Herrmann, M. Karweit, G. Drazer, Phys. Rev. E 79 (2009) 061404. 

[35] J. Frechette, G. Drazer, J. Fluid Mech. 627 (2009) 379-401. 

[36] M. Balvin, E. Sohn, T. Iracki, G. Drazer, J. Frechette, Phys. Rev. Lett. 
103 (2009) 078301. 

26 



[37] H. Wu, W. N. Everett, S. G. Anekal, M. A. Bevan, Langmuir 22 (2006) 
6826-6836. 



[38 
[39 

[40 

[41 
[42 
[43 

[44 

[45 
[46 

[47 
[48 
[49 



P. Bahukudumbi, M. Bevan, J. Chem. Phys. 126 (2007) 244702. 

H. Brenner, D. A. Edwards, Macrotransport Processes, Butterworth- 
Heinemann, Boston, 1993. 

H. Risken, The Fokker-Planck Equation: Methods of Solutions and Ap- 
phcations. Springer, Berhn, 1989. 

M. H. Jacobs, Diffusion Processes, Springer, New York, 1967. 

R. Zwanzig, J. Phys. Chem. 96 (1992) 3926-3930. 

P. S. Burada, G. Schmid, D. Guerrera, J. M. Rubi, P. Hanggi, Phys. 
Rev. E 75 (2007) 051111. 

N. Laachi, M. Kenward, E. Yariv, K. D. Dorfman, Europhys. Lett. 80 
(2007) 50009. 

X. Wang, G. Drazer, Phys. Fluids 21 (2009) 102002. 

W. B. Russel, D. A. Saville, W. R. Schowalter, Colloidal dispersions, 
Cambridge University Press, USA, NY, 1989. 

D. L. Ermak, J. A. McCammon, J. Chem. Phys. 69 (1978) 1352-1360. 

Y. Pawar, J. L. Anderson, Ind. Eng. Chem. Res. 32 (1993) 743-746. 

M. E. Staben, A. Z. Zinchenko, R. H. Davis, Phys. Fluids 15 (2003) 
1711-1733. 



27 



